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We report on high resolution measurements of resonances in the spectrum of coherent synchrotron 
radiation (CSR) at the Canadian Light Source (CLS). The resonances permeate the spectrum at 
wavenumber intervals of 0.074 cm -1 , and are highly stable under changes in the machine setup 
(energy, bucket filling pattern, CSR in bursting or continuous mode). Analogous resonances were 
predicted long ago in an idealized theory as eigenmodes of a smooth toroidal vacuum chamber 
driven by a bunched beam moving on a circular orbit. A corollary of peaks in the spectrum is 
the presence of pulses in the wakefield of the bunch at well defined spatial intervals. Through 
experiments and further calculations we elucidate the resonance and wakefield mechanisms in the 
CLS vacuum chamber, which has a fluted form much different from a smooth torus. The wakefield 
is observed directly in the 30-110 GHz range by RF diodes, and indirectly by an interferometer in 
the THz range. The wake pulse sequence found by diodes is less regular than in the toroidal model, 
and depends on the point of observation, but is accounted for in a simulation of fields in the fluted 
chamber. Attention is paid to polarization of the observed fields, and possible coherence of fields 
produced in adjacent bending magnets. Low frequency wakefield production appears to be mainly 
local in a single bend, but multi-bend effects cannot be excluded entirely, and could play a role 
in high frequency resonances. New simulation techniques have been developed, which should be 
invaluable in further work. 
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INTRODUCTION 

A relativistic particle bunch moving along a curved tra¬ 
jectory within a conducting metallic chamber is accom¬ 
panied by an electromagnetic wake generated by the in¬ 
teraction between the bunch and the chamber. When the 
bunch or a substructure within it is sufficiently short, the 
particles can act coherently to produce the wake, and the 
emitted power varies as the square of the current. This 
implies an enormous intensity compared to the usual in¬ 
coherent SR at the same frequency. The character of the 
wake is strongly influenced by the chamber. For exam¬ 
ple, the wake of a bunch moving along a circle of radius 
R between parallel plates of separation h is weakened 
due to the shielding action of the induced image charges 
0, and as a result the CSR is exponentially suppressed 
at wavelengths greater than Ao ~ 2 h(h/R) 1 / 2 , called the 
CSR shielding cutoff. 

The situation changes dramatically when conducting 
side walls are added to form a torus [2| . Now the spectral 
distribution of CSR is predicted to have a series of narrow 
peaks associated with excitation of resonant modes of the 
torus, and the shielding cutoff can be at a longer wave¬ 


length, depending on the radial beam-to-wall distance. 
A transfer of power into these modes occurs when the 
phase velocity of the modes matches the particle veloc¬ 
ity. There is a close analogy to whispering gallery modes 
in acoustics, which are described by similar mathemat¬ 
ics 0- 

The interaction between the field and the bunch can 
be characterized by an impedance, which we define in 
a more general sense than usual. Suppose that the 
purely longitudinal current density can be factored as 
I(s — ct)f(x,y) where ( s,x,y) are the usual accelera¬ 
tor coordinates (arc length on a reference trajectory and 
transverse displacements). A general representation of 
the i— th component of the full field is Ei(s,x,y,t ) = 
— f exp(2wik(s — ct))Zi{s,x,y,k)I{k)dk , where Zi is 
the *-th component of the impedance. We define the 
wavenumber k = 1/A without a factor of 2-7T. The 
position-dependent wakefield, the object of the present 
study, is the full field regarded as a function of the dis¬ 
tance from the bunch, namely £i(z, s , x, y ) = Ei(s , x, y, t) 
with z = s — ct. We study the time dependence of 
£ x (s — ct, r), and confirm that it depends on the position 
r = {s,x,y). 
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Peaks in the impedance as a function of k correspond to 
resonances. It is expected that resonances in Z x will cor¬ 
relate with those in Z s , as happens in the analytic theory 
[§j. Structures complementary to resonances appear in 
the wakefield, through the Fourier transform. Consider, 
for example, an idealized torus impedance consisting of a 
train of resonances uniformly spaced by A / in frequency. 
The corresponding time domain wakefield consists of a 
train of localized pulses, equally spaced by the distance 
Az given by the approximation 
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where D is the distance from the orbit to the outer wall 
of the torus. This periodicity is consistent with detailed 
wakefield calculations by Agoh [dj, Stupakov and Kotel- 
nikov Q, and Zhou et al. p], for radiation in rectangu¬ 
lar chambers of various cross sections, usually in a sin¬ 
gle bend. Particularly noteworthy is the highly localized 
nature of the wakefield pulses predicted by these calcu¬ 
lations. 

In this report we present our time domain measure¬ 
ments of the wakefields generated within the vacuum 
chamber in the bending magnets (R = 7.143 m, bend 
angle 0 = 15°, h = 32 mm) in the Canadian Light Source 
(CLS) storage ring. The wakefield following the bunch is 
indeed found to consist of a sequence of localized pulses, 
but the pattern is less regular than in the toroidal model, 
and depends on the point of observation in the vacuum 
chamber. At the higher frequencies (6 - 15 cm -1 , 180 - 
450 GHz) the wakefield was characterized indirectly using 
a high resolution Michelson interferometer. At lower fre¬ 
quencies (1.1 - 3.7 cm -1 , 33 - 110 GHz) backward prop¬ 
agating fields were observed using a series of RF diodes 
mounted on a forward-viewing port of the vacuum cham¬ 
ber, referred to as backward detectors. Note that by their 
nature our detectors respond only to transverse compo¬ 
nents of the fields. 


INTERFEROMETER MEASUREMENTS 

The Michelson interferometer measures the interfer¬ 
ence between a direct and a delayed far infrared (FIR) 
beam, and can observe peaks of intensity when the 
path length difference is a multiple of the spacing be¬ 
tween pulses in the wakefield. The interference pat¬ 
tern is given by 2 A 1 A 2 (E x (t + r)E x (t)), where Ai and 
A 2 are amplitudes including propagation losses of the 
two beams in the instrument, AiE x is the horizon¬ 
tal electric field, and r is the time delay, with (•) 
representing the average over t Q. Expressing the 
field in terms of the horizontal impedance and av¬ 
eraging over a long time we see that the intensity 
is proportional to f exp(2irikcT)\Z x (r, k)I(k)\ 2 dk, an 


even function of r. The frequency domain current 
I(k) may contain contributions from several bunches. 
By contrast the wakefield of interest is proportional 
to f exp(2TrikcT)Z x (r, k)I(k)dk, but the two quantities 
should both show peaks in r at roughly the same points. 

Radiation for the interferometer measurements was 
produced by operating the ring at 2.9 GeV, with the mo¬ 
mentum compaction adjusted to give a bunch length of a 
few millimeters. The CSR was produced in the bursting 
mode by running at a current of 12 rnA distributed uni¬ 
formly over 3 consecutive bunches. Coherence of the ra¬ 
diation was confirmed by noting a quadratic dependence 
on the bunch current. 

Spectra were acquired on the FIR beamline using a 
Bruker 125 HR interferometer with a 75 /im Mylar beam 
splitter and an Infrared Laboratories Si bolometer as de¬ 
tector with a gain at 200X. We used a 12.5 mm aperture, 
a scanner velocity of approximately 0.051 m/s (80 kHz as 
measured using a HeNe laser reference), and a resolution 
of 0.020 cm -1 (defined as 0.9 / path-length-difference by 
Bruker convention). The spectra fall mostly in the fre¬ 
quency range 6-15 cm -1 or 0.18-0.45 THz. The intensity 
versus path length difference is shown in Fig |T] The sig¬ 
nal comes through the FIR beamline from the Ml mirror 
at the end of the bend. Of special note is the strong 
interference pattern at 13.5 cm and the weaker patterns 
near 27 cm and 41 cm, which we associate with trailing 
wakes as discussed in the Introduction. We have verified 
that the interference patterns disappear with incoherent 
synchrotron radiation, which helps to rule out instrumen¬ 
tation anomalies as a factor in the measurements. 



FIG. 1: Interferogram as a function of path length difference 


The polarization of the signal transmitted by the FIR 
beamline was analyzed using a wire grid polarizer with 
a wire spacing of 10 /im, placed at the entrance to the 
interferometer. The vertical polarization (as it would ap¬ 
pear at the Ml mirror) was found to be negligible, less 
than 0.5% of horizontal. 

Turning now to the Fourier transform of the interfer¬ 
ogram, we find the result of Fig. [2] The fine structure 
is spaced by A k ~ 0.074 cm -1 , the reciprocal of the 
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13.5 cm spacing of Fig. [T] We have accumulated simi¬ 
lar data over a period of years, and of particular note 
is the stability of the spectral pattern of Fig. [2], includ¬ 
ing phase, under wide variations in the electron beam 
parameters such as energy (1.5 or 2.9 GeV), fill pattern 
(single bunch or 210 bunches), and CSR mode (burst¬ 
ing or continuous). Also, large changes in the beamline 
optics were made over the years. A graph showing this 
invariance appears in Ref. [3]. This stability argues in fa¬ 
vor of the vacuum chamber as being the dominant factor 
in shaping the spectrum. 



FIG. 2: Fourier transform of the interferogram 


RF DIODE MEASUREMENTS 

The vacuum chambers in the bending magnets are 
equipped with quartz vacuum windows on 50 cm long 
tubular extensions mounted near the chamber entrance 
(Fig. [3]). These align with the axis of the downstream 
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FIG. 3: (color online). Fluted vacuum chamber at the FIR 
dipole with bending radius R = 7.143 m and deflection angle 
8 = 15°. The maximum excursion of the outer wall from 
the beam is 33 cm. The diode detector is at the end of the 
horizontal pipe at the left of the figure. 

straight section. RF Schottky diodes (unbiased) mounted 
outside the windows have proved to be highly sensitive 
to backward-propagating fields in the 30-110 GHz range. 
These fields are correlated in time with the bunch and 
appear to be caused by back-scattering of waves by hard¬ 
ware near the chamber exit. This hardware consists of 


two parts: a bar-like copper photon absorber which ex¬ 
tends about 20 cm into the chamber, at a right angle to 
the outer wall, and the Ml mirror assembly which di¬ 
rects the FIR radiation upward into the FIR beamlinc. 
Three diodes were used spanning 33-50 GHz, 50-75 GHz 
(with a 50-75 GHz band-pass filter), and 75-110 GHz 
[8:j. The detector waveguides determine the bandwidths 
and polarization axes. To the best of our knowledge, the 
first application of Schottky diodes to the measurement 
of freely propagating wakefields is reported in Ref Q. 

CSR for the diode work was produced at 2.9 GeV in the 
continuous mode, with a single bunch of length 2 mm (as 
deduced from the synchrotron frequency) carrying 1 mA. 
The diodes are highly sensitive to polarization, and we 
have used this feature to separate the electric fields polar¬ 
ized parallel and perpendicular to the orbital plane, com¬ 
plementing the interferometer study. The signals from 
the 50 - 75 GHz diode for these two polarizations are 
presented in Fig. [4] as curves 1 and 2 respectively. Our 
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FIG. 4: (color online). RF diode measurements in the 
time domain (oscilloscope traces) with a 50-75 GHz detector. 
Diode mounting and polarization: 1-backward horizontal; 2- 
backward vertical; 3-forward horizontal (with adjustment of 
time base). For clarity the curves have been separated verti¬ 
cally. 

interpretation of the figures is as follows. The horizontal 
data (Curve 1) show several downward peaks which we 
partition into two groups consisting of the first pair (A- 
B) and the cluster (C-G). The two groups correspond to 
reflections from the photon absorber bar and the Ml mir¬ 
ror assembly, respectively, the latter being downstream 
of the former. Peaks A and C are interpreted as the 
prompt signal from the bunch as reflected from the pho¬ 
ton absorber and the mirror, respectively. Accordingly, 
the A to C distance is 46 cm (1.53 ns), which is twice the 
physical separation of the two structures (measured cen¬ 
ter to center). This interpretation is supported by a test 
measurement in a bending magnet which had no mirror 
assembly and only peaks A and B survived. 
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Peak B is understood as a wake pulse from the pho¬ 
ton absorber following A, whereas D, E, F, G are wake 
pulses from the mirror following C. One might ask why 
there are no further wake pulses from the photon ab¬ 
sorber following B, but that is exactly the behavior seen 
in a simulation reported presently. Finally, the remaining 
structures in Curve 1 immediately preceding the peak C 
correlate with absorbers protecting the mirror assembly 
and are of no further interest here. 

The results for vertical polarization, in Curve 2, show 
little intensity prior to signals reflected from the mirror 
assembly. This makes sense, because the mirror assem¬ 
bly is the first structure with strong asymmetry about 
the midplane y = 0. If electromagnetic boundary con¬ 
ditions (and the charge distribution) were symmetrical, 
the predicted vertical polarization would be zero. 

The RF diode measurements have been repeated at the 
exit of the FIR beamline, directly in front of the interfer¬ 
ometer, and results for horizontal polarization are shown 
in Fig. [f]as Curve 3. With the time base adjusted to line 
up the prompt signal (the large peak) with peak C, the 
subsequent patterns observed at forward and backward 
detectors line up rather well. The finer oscillations pre¬ 
ceding the big peak in Curve 3 are a manifestation of the 
approximately 3.5 GHz bandwidth of the oscilloscope. 

Although the diode and interferometer measurements 
are not strictly comparable, operating as they do in dif¬ 
ferent ranges of frequency, the spacing of wakefield pulses 
as transmitted from the Ml mirror appears to be roughly 
the same in the two measurements. For instance, the 
spacing in Curve 3 of Fig.[4]is about 15 cm, versus 13.5 cm 
in Fig. |H 

We turn to the question of wake propagation between 
bending magnets. Are wakes produced locally or is there 
cooperation between adjacent dipoles? A cooperative ef¬ 
fect was computed in @ and [fl] for a chamber of constant 
cross section, but it is parameter dependent and the prob¬ 
lem is complicated by our fluted geometry. For an exper¬ 
imental test a different bending magnet was used, one 
5.4 m downstream from an adjustable vertical scraper. 
Otherwise the setup is identical to the FIR dipole, and 
the backward diode signal is identical to that of Curve 1 
of Fig. [f] The sum of the intensities of leading peaks 
A and B was measured as a function of vertical aper¬ 
ture over the range 30-2 mm. If a large part of the 
A + B intensity came from a bend or bends before the 
scraper, one might expect the intensity to fall markedly 
when the scraper aperture is at the minimum, at least 
if the propagating mode were broadly distributed over 
the open aperture as in the lowest vertical mode in a 
waveguide. Instead we found that the A + B intensity 
increased by about 10%, while the form of the peaks was 
more or less invariant. This increase is presumably due to 
enhanced CSR in the local bend due to a perturbation of 
the bunch form by the scraper impedance. This indicates 
that the A + B wake is produced locally to a consider¬ 


able extent, since in the straight between the scraper and 
the local bend there is nothing to directly amplify the 
radiation. On the other hand, if only a portion of the 
final wake was produced upstream and was attenuated 
by the scraper, the attenuation could be masked by the 
enhanced local CSR. Thus we fail to exclude entirely a 
partial production of the wakefield in upstream magnets. 

SIMULATING THE FIELD PATTERN 

To simulate fields in the fluted chamber we made a 
time-domain integration of the six Maxwell curl equa¬ 
tions, with boundary conditions, by the Discontinuous 
Galerkin Method [lOj. A Fourier expansion in y leaves 
two spatial coordinates (s, x). The receding chamber wall 
is modeled by a flat conducting surface, connected at its 
end to the photon absorber represented as a rectangle. 
After the absorber, the wall runs parallel to the beam at 
x = 33 cm and then almost perpendicular to the beam 
to meet the mirror assembly, which is rendered as a flat 
surface, tipped by 7.5° to aim at the center of the bend. 
A Gaussian driving bunch with cr = 2 mm excites fields 
in a single pass, the initial condition being the steady 
state fields in a straight wave guide. A straight exit 
channel after the mirror is long enough to prevent re¬ 
flections from the end of the channel reaching the bend 
region. We used 61195 finite elements of 6th order, 28 
nodes per element, for a maximal node spacing of about 
1 mm. The time step was cdt = 0.16 mm. The pipe to 
the diode is not modeled, but the time dependence of the 
fields is recorded at the pipe entrance, here denoted as 
the “backward port”. The diode voltage for horizontal 
polarization is proportional to E%, as conditioned by the 
detector response. We account for the latter in the fre¬ 
quency domain by the low pass filter 1/(1 + if/fo), with 
/o = 1 GHz. 



FIG. 5: (color online). Simulated E% at backward port vs. ct, 
after a low pass filter to account for detector response. The 
origin of time t is when the bunch is 5 cm before the entrance 
to the bend. Only the lowest mode in y is included. 

The calculated detector response at the backward port 
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shown in Fig J5] displays signals reflected from the photon 
absorber and the mirror, just as in the picture inferred 
from experiment before the simulation. Field contour 
plots show sharply defined wave fronts in the reflected 
fields. The red lines in Figj5] mark the positions of the 
experimental peaks in FigH] (relative to the first theo¬ 
retical peak). In view of our rudimentary modeling of 
the reflecting structures, which actually have significant 
three-dimensional aspects, the resemblance to the experi¬ 
mental wake pulse positions seems quite satisfactory. The 
absence of a third wake pulse in the wave from the pho¬ 
ton absorber is striking, and emphasizes the dependence 
of the wake pattern on the point of observation. The cal¬ 
culated pulse D is much too strong and a bit too early, 
perhaps due to over-simplified modeling of the mirror 
structure. 

The unfiltered signal shows significant structures be¬ 
tween peaks A and B, and between C and D. These might 
be related to the pulse at 7 cm in the interferogram. 

The code was validated by checking boundary condi¬ 
tions and by comparison with a well tested paraxial code, 
which works at points not too close to the end of the 
bend, before any backward waves are excited. A simula¬ 
tion of the interferometer data, lying at higher frequen¬ 
cies, will require higher resolution or possibly a frequency 
domain method. 
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